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The origin of high velocities of pulsars is studied by considering the spin-flip conversion of neutrinos 
propagating in a gravitational field of a protoneutron star. For a rotating gravitational source (such 
as pulsars) with angular velocity u>, one finds that the spin connections (entering in the Dirac 
equation written in curved space time) induce an additional contribution to neutrino energy which 
is proportional to uj p, with p the neutrino momentum. Such a coupling (spin-gravity coupling) can 
be responsible of pulsar kicks being the asymmetry of the neutrino emission generated by the relative 
orientation of the neutrino momentum p with respect to the angular velocity ui. As a consequence, 
the mechanism suggests that the motion of pulsars is correlated to their angular velocity us. In 
this work we consider neutrinos propagating orthogonally to the magnetic field. The fractional 
asymmetry turns out to be independent on the magnetic field of the nascent protostar, and is only 
related to the angular velocity (Ap/p ~ u>). As in the usual approaches, spin flip conversion is 
generated via the coupling of the neutrino magnetic momentum with the magnetic field. For our 
estimations, we use the large non-standard neutrino magnetic momentum provided by astrophysical 
and cosmological constraints, [i v ~ I0 _11 /is. 

The connection with recent observations and statistical analysis is also discussed. 



I. INTRODUCTION 

A discussed and unsolved issue of the modern astro- 
physics is the origin of the pulsar velocity, i.e. the high 
proper velocities of pulsars as compared with the sur- 
rounding stars. Their three-dimensional galactic speed 
runs, in fact, from 450 ± 90Km/sec up to a maximum of 
about lOOOKm/scc lj. This peculiarity immediately has 
suggested that nascent pulsars receive an impulse (kick) 
during their formation. The gravitational binding energy 
(3 x I0 53 erg) is, after the supernova collapse of a massive 
star, carried out by outgoing neutrinos (about 99%). An 
anisotropy of ~ 1% of the momenta distribution of the 
outgoing neutrinos would then suffice to account for the 
neutron star recoil of 300Km/sec. 

An interesting mechanism to generate the pulsar ve- 
locity has been recently proposed by Kusenko and Segre 
(KS) 0- It involves the physics of neutrino oscillation in 
presence of an intense magnetic field. Let us recall the 
basic idea. The neutrinosphere is defined as the surface 
from which neutrinos may escape from the protostars. In 
particular, electron neutrinos v e are emitted from a sur- 
face which is located at a distance from the center greater 
than the surfaces corresponding to muon/tau (^,r) neu- 
trinospehres. Under suitable conditions, a resonant os- 
cillation v e — > i/^t can occur between the v e and f MiT 
neutrinospheres. Neutrinos generate via oscillations 
can escape from the protostar being outside of their neu- 
trinosphere, with the ensuing that the "surface of the res- 
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onance" acts as an "effective f MjT -neutrinosphere" . The 
presence of a magnetic field may distort the effective sur- 
face of resonance and the energy flux turns out to be 
generated anisotropically. In the KS mechanism Q, the 
responsible for the neutrino emission anisotropy is the 
polarization of the medium due to the magnetic field B. 
The usual MSW resonance conditions turn out to be, 
in fact, modified by the term (first derived by D'Olivo, 
Nieves and Pal) 0] 

eG F /3n e \ 1/3 „ „ 

where p = p/p, p is the neutrino momentum, e is the 
electric charge, Gf is the Fermi constant, and n e is the 
electron density. The KS mechanism has been also stud- 
ied for active-sterile neutrino oscillations (sterile neutri- 
nos may have a small mixing angle with the ordinary 
neutrinos) |4j, |5j . Papers dealingwith the origin of pul- 
sar kicks can be found in Jl 0, HJ1M IH 111 13 III 

III HI III El El El El H illlU 

In this paper is discussed the possibility that spin fla- 
vor conversion of neutrinos propagating in a gravitational 
field generated by a rotating source may generate the ob- 
served pulsar kicks [24[ (the role of the angular velocity 
on pulsar kicks has been also studied, although in a differ- 
ent context, by Mosquera Cuesta 18]). The gravitational 
field affects, as we will see, the resonance conversion of 
left-handed neutrinos into right handed neutrinos, the 
latter being sterile can escape from the neutrinosphere, 
hence from the protoneutron star. In particular, the 
modification to the resonance condition is induced by 
spin connections which enter in the Dirac equation in 
curved spacetimes. They give rise to a coupling term 



2 



~ U3 ■ p, where u) is the angular velocity of the gravi- 
tational source. The relative orientation of neutrino mo- 
menta with respect to the angular velocity determines an 
asymmetry in the neutrino emission, hence may generate 
pulsar kicks. 

The paper is organized as follows. In Sect. II we re- 
view the Dirac equation in curved space-times. Here we 
see that, owing to the breakdown of the spherical symme- 
try generated by the angular velocity of the gravitational 
source, spin connections are non vanishing and are pro- 
portional to the chiral operator 7 s . Sect. Ill and Sect. IV 
are devoted, respectively, to briefly recall the main fea- 
tures of matter induced effective potential and the role 
and intensity of magnetic fields in astrophysical systems. 
In Sect. V the fractional asymmetry is computed. Here 
we also discuss the resonance and adiabatic conditions, 
as well as the spin flip probability that left-handed neu- 
trinos transform in right-handed neutrinos. Conclusions 
are drawn in Sect. VI. 



where F A1 (x) are the spin connections 



(semicolon represents the covariant derivative). Using 
the relation 

dr b ci rthb c r,hc b o - ddbc 5 
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where r} & t is the metric of flat spacetime, and after some 
manipulations, the spin connections can be cast in the 
form 26] 

yiV = 7 a <rv = 7 a < [iAap [-(-g)- 1/2 7 5 ] } , 



where 
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II. DIRAC EQUATION IN CURVED SPACE 
TIME 

The phase of neutrinos propagating in a curved back- 
ground is generalized as [2€ 



l^/(A)) 
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(1) 



where / is the flavor index and j the mass index. Utj are 
the matrix elements transforming flavor and mass bases 



U = 



cos tf sin a 
— sin 6 cos 



(2) 



(in what follows we shall consider the neutrino mixing 
for two flavors). Besides, 

P ■ Pnull = P^Pnull 1 

where is the four-momentum operator generating 
space-time translation of the eigenstates and 



fi 

Pnull — 



is the tangent vector to the neutrino worldline x^, pa- 
rameterized by A. From here we follow the Cardall and 
Fuller paper |26j . The covariant Dirac equation in curved 
space-time is (in natural units) 28] 

[i^{x)D ll - m]ip = 0, 

where the matrices j^(x) are related to the usual Dirac 
matrices 7° by means of the vierbein fields e^(x), i.e. 

7^)=e£(*h a . 

The Greek (Latin with hat) indices refer to curved (flat) 
space-time. is defined as 



D„ 



g = det(.g M „) . 

The above procedure allows to separate out the chirality 
operator 7 s . This shows that 7 M r M acts differently on 
left- and right-handed neutrino states. In fact, by writing 

Pp - Pl , 
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where 



P, 



L,R 



1T7 5 



are the projection operators, one sees that neutrinos with 
left- and right-handed helicity acquire a different gravita- 
tional contribution. In the case of neutrino oscillations, 
one can add, without physical consequences, a term pro- 
portional to the identity matrix (~ AguI), so that 7 s 
can be replaced by the left-handed projection operator 
Vl = (1 — 7 5 )/2- As a consequence, the spin-gravity 
coupling is pushed in the left-handed sector of the effec- 
tive Hamiltonian (see Eqs. ©-©), and no contributions 
appear in the right handed sector. 

The equation of evolution of neutrinos has the form 
(/ = e, f' = n,T) 







( VfL 


Vf'L 


= H 


v f'L 


VfR 




1 VfR 
\ Vf'R 
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dX 

where the diagonal terms of the Hamiltonian TL 

diagW = U~ l MU + Q G {x) V L , 

are written in terms of the mixing matrix U and mass 
matrix M 



M = 



m\ 
ml 
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and fl G is defined as 



E 



•mi and m 2 are the mass eigenstates, 9 is the vacuum 
mixing angle, and = {E, p), being E the energy mea- 
sured in the local frame. 

The inclusion of matter induced effective potential and 
magnetic terms, which appear in the off-diagonal terms 
of the Hamiltonian H, will be discussed later. 



A. The Geometry of a Rotating Mass Source 

For geometries with a spherical symmetry, such as the 
Schwarzschild or Reissner-Nordstrom space-times, A G 
vanish. Nevertheless, for rotating gravitational sources 
(Lense-Thirring geometry), whose line element is (in 
weak field approximation) 

ds 2 = (1 - (p)(dt) 2 - (1 + 4>){d-x) 2 -2h-d*dt, (3) 

with 



(x, y, z) , r — y/ x 2 + y 2 + z 2 



, 2GM , / 4GM R 2 \ 
#»•) = — > h =[ 5r 3 j^Ax, 

u> is the angular velocity of the gravitational mass M, 
and R its radius, A G acquire a non vanishing component 
related to the off-diagonal terms of the metric tensor, 
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where 



uj = u> — 
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One can show the angular velocity induces a drift velocity 
of neutrinos p9|. 

The vierbein fields used for computing the spin con- 



nection A G are 



As the angular velocity goes to zero, the spherical sym- 
metry is recovered and the spin connections vanish, as 
immediately follows from the above expressions. 

It is worth note that the term 7 A T M can be rewritten 
in the form {y/—g ~ 1 + 2(f)) 



17% ~ i^{'iA G ^) = /(r) 7 °o;" i I]' i 



where 



and 
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Here a' 1 , i = 1,2,3 are the Pauli matrices. 

The term ~ /(r)o;' • S is the well known 
gravitomagnetic-spin coupling (see for example |3Cj). 



III. MATTER INDUCED EFFECTIVE 
POTENTIAL 

Neutrinos inside their neutrinospheres are trapped ow- 
ing to weak interactions with the background matter, 
which lead to the potential energy 



Vf 



3.8 



P 



10 14 gr cm 



zioVfir, t)eV, 



where p = m n n e is the matter density, m n is the nucleon 
mass, f = e,fi, t, and 



Ve 



n-i/3, 



In these expressions, Y e is the electron fraction. In the 
present analysis we shall envisage those neutrinos for 
which matter induced effective potential for (electron) 
left-handed neutrinos such that V Ve <C 1. This follows 
in the regions where the electron fraction Y e assumes the 
value sa 1/3 (v f <g 1) [H l3lL These regions are 

located at r ~ 15 km (see the paper by Nunokawa, Pel- 
toniemi, Rossi and Valle |3*^1. 



IV. MAGNETIC FIELD INSIDE A PROTOSTAR 



e ol = -h i , e i6 = 0. 

The non vanishing A G is an indication of a preferred 
direction related to the angular velocity of the source. 
n G (x) can be rewritten as 



E, 



5r 3 Ei 



(4) 
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Mq VIOkmy \ r J 10 4 Hz 



Electroweak interactions of neutrinos with matter 
background play a central role on the neutrino emission 
during the core collapse of supernovae, and in general, 
on star cooling mechanism of magnetized medium. The 
neutrino energy spectrum in presence of strong magnetic 
fields is modified, depending on their flavors. 

Protoneutron stars possess strong magnetic fields 
whose strength is > 10 12 — 10 14 G. Such strong fields 
can be also found near the surface of supernovae [34| . 
neutron stars |35j, and magnetostars 36]. A fundamen- 
tal feature of large magnetic fields in such astrophysical 



4 



systems is related to their effects on neutrinos, which 
through charged and neutral current interactions, mod- 
ify the dispersion relations of neutrinos. 

In computing the effects of strong magnetic fields on 
neutrinos propagation, the following condition holds 

m 2 e < eB < Mw, 

where m e is the electron mass and Myy is the W-boson 
mass. We define the fields 

2 

B c = — ~ 4.4 x 10 13 G ~ 8.62 x 10 n eV 2 . 

e 

It "separates" the regimes of weak field B -C B Cl and of 
strong field B > B c . 

In a series of recent papers j33, HU , a detailed analysis 
of the neutrinos physics in strong magnetic fields has been 
carried out. In what follows we shall consider those neu- 
trinos which propagate orthogonally to the magnetic field 
of the protostar, so that B • p = Bcosa = (a = tt/2). 
The only contribution to the neutrino energy comes from 
the magnetic momentum of neutrinos. 



where F~; is the electro- magnetic field tensor, and 

ab ° 5 

hi 1 r a _,Sl 

a = jil ,7 J • 

In this paper, we shall consider the neutrino magnetic 
momentum 

Hf ~ (10~ 12 - 10- n )/i B , 

as provided by astrophysical and cosmological constraints 
|40|. Such a value on the neutrino magnetic momentum 
is great as compared with ones of the Standard Model 
prediction fiff ~ 10~ 19 LtB, provided that the neutrino 
mass is m„ ~ few eV |41| . 

The expression for the neutrino magnetic energy is 

Uff>B ~ 2.5 x iq-6 Vff' _g_ eV , 
W 10- n /iB B c 

V. EQUATION OF EVOLUTION OF 
NEUTRINOS AND THE ASYMMETRIC 
NEUTRINO EMISSION 



A. Neutrino Magnetic Momentum 

Since neutrinos are uncharged particles, they do not 
interact directly with photons (magnetic fields). The 
typical coupling of the electromagnetic field with the 
fermionic current is absent for neutrinos. This is essen- 
tially due to the fact that the Standard Model is built 
up assuming that neutrinos are massless, hence only left- 
handed neutrinos appear in the theory. The experimental 
evidence of neutrino oscillations is an index that neutri- 
nos are massive particles. In such a circumstance, we 
are in the context of a theory which goes beyond the 
Standard Model since right-handed projection of neutri- 
nos has to be included in the fermion sector. As a con- 
sequence, an anomalous magnetic momentum emerges 
through quantum corrections (one-loop diagram whose 
internal lines are charged lepton and W-Z boson propa- 
gators). Thus, even though neutrinos have no charged, 
they posses a magnetic momentum, which induces the 
interaction with photons 



10 



-19 



eV 



where m v is the neutrino mass, and 

e 

fl>B 



2m f 



is the Bohr magneton. 

Neutrinos interacting with the magnetic field of the 
protoneutron star acquire, hence, an energy through the 
interaction 391 



Taking into account for the gravitational and magnetic 
contributions, the equation of evolution describing the 
conversion between two neutrino flavors / = e and /' = 
fx, t reads 01 



dX 



( V fL 

VfR 
\ VfR 




(5) 



where, in the chiral base, the matrix Ti is the effective 
Hamiltonian defined as 



H 



H L H fJ , 
Hff> Hr 



5c 2 



Ss 2 V U f, + + Sc 2 



H 



—Sc2 Ss 2 
Ss 2 Sc 2 



(iff, fi f , f , 



(0) 



(7) 



(8) 



6 = 



AE 



km 1 



m 9 



Ant = lpflv(T &b F &i) llj , 



c 2 = cos 29 , s 2 = sin 29 , 
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B± = Bs'ma = B is the component of the magnetic 
field orthogonal to the neutrino momentum, and (3 is the 
angle between the neutrino momentum and the angular 
velocity. 

Qg is diagonal in spin space, so that it cannot induce 
spin-flips[54(. Its relevance comes from the fact that it 
modifies the resonance conditions of spin-flips 

f/L V Ve +Sl G (f)-25c 2 = 0, (9) 
Vf'L -> VfR V Vf , + n G (f) + 2Sc 2 = . (10) 

f is the radial distance where the resonance occurs. The 
resonances (Q and I|1U|) do not occur simultaneously; in 
what follows we shall consider the transition (JHJ). 

Besides, we shall use the best fit for solar neutrinos 

El 



Am 



Sun 



(1(T 5 l(T 4 )eV 2 



sin 2 29 S un > 0.8 . (11) 



fta in I® (and l|lUfl) distorts the surface of resonance 
due to the relative orientation of the neutrino momentum 
with respect to the angular velocity. As a consequence, 
the outgoing energy flux F results modified. The neu- 
trino momentum asymmetry is defined as |ll| 



Ap _ 1 fo F • uda 
p 3 JqF- hda 



(12) 



where the factor 1/3 accounts for the conversion of neu- 
trinos Vfi, into VfR (f 7^ /'), da is the element of area 
on the distorted surface, u) and n are the unit vectors 
parallel to the angular velocity and orthogonal to da, 
respectively. To compute the fractional asymmetry one 
then should specify the protostar model |44|. including 
into the hydrodynamical equations the rotational effects 
due to gravitational sources j^. This task goes beyond 
the aim of this work. 

Here we give an estimation of the fractional asymme- 
try following the calculations as in KS Q ■ The surface of 
resonance is located to radial distance r((3) = f + gcos (3 
(cos/3 = p • u)', with cl>' = lj'/lj'). Inserted in the reso- 
nance equation, one has 



f(r) 



V' 



m 

V 



where we used the resonance condition 



2Sc 2 = V Ve (r), 



(13) 



(14) 



and 0' G < Vl t 



as evaluated at f. The prime means 
derivative with respect to r. In such a calculations we 
assumed that the matter density profile is described by 

M 



P(r) = pa 



where po = 2 10 14 gr/cm 3 and p = 0(1). The fractional 
asymmetry reads 



Ap 

P 



2g dT 
3T~d7 



2/(r) (dV Vc 



= 82.5 



3 T V dT 
MeV 3 /(f) (dn^ ' 
Ty e ^V" \lT 



(15) 



Since the density number n e is related to the temperature 
T via the Fermi distribution 



n P = 2 



d 3 p 



1 



(2tt) 3 e(P"^)/ T + 1 ' 



(16) 



where /x is the chemical potential, one obtains dn e /dT ~ 
2r]T 2 /3, where 77 = fx/T ~ 5 (see the paper by Qian Q). 
Eq. (|15|l then becomes 



Ap 

P 



0.01 



4.5 10" 



/3MeVV 1 /(f) 



(17) 



V T J r)y e eV 

From the resonance condition one obtains 

Am 2 cos 2(9-7.6 10 4 y e eV 2 , 

as the density p ~ 10 11 gr/cm 3 at f ~ 15km. The value 
falls into the best fit (|llf> provided y e w 10~ 9 . One can 
then see that for typical values of pulsars 



i?>15km, M ~ M, 



■ 



T - 3MeV , 



Eq. I|17|) implies that the fractional asymmetry — 1% is 
recovered as 

uj - few 10 2 Hz , 
which is, indeed, the typical angular velocity of pulsars. 

A. Adiabatic Conditions 

Besides the level crossing it must be also adiabatic, 
i.e. the corresponding (adiabatic) parameter 7, which 
quantifies the magnitude of the off-diagonal elements 
with respect to the diagonal ones of © in the instan- 
taneous eigenstates, must satisfy the condition 7(f) 3> 1. 

The precession length L coming from the effective 
Hamiltonian JSJ is given by ^7j 



L = 



2tt 



^(2pi fr B ± )i + (K e + n G - 2Sc 2 y 

At the resonance, it reads 



L r es — L(f) — 



7T IO^Vb B c 
— — ~ 87T m . 

fj,ff>B± [i f f, B± 



The resonant spin flavor precession has the following 
width 

Sr = 2 AX , 



6 



where 



VI. CONCLUSIONS 



A 



P'^Y 1 V Ve (r) 



p(r) 



V'(r)' 



assuming 



and 



y' f (t,t)=0, (y e '«n' e /n e ) 



A = 2- 



The spin flavor conversion is adiabatic provided 

Sr > L res , 

which is equivalent to 
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This condition can be immediately verified using the 
above profile for p(r), and (iff ~ 10 _1 Vb and the typ- 
ical values for the magnetic fields B ~ (10 12 — 10 14 )G. 



B. Spin Flip Probability 

The conversion probability P UfL — n/ /R that the neu- 
trino left flips into neutrino right is 

Pu SL ^u f , R = \ - Q - PLZ^j COS 29% COS 29 f , 

where the Landau-Zener probability is given by 

Plz = e-^' 2 , 
and the effective mixing angle 9 is defined as 0, 

2p sr B 



tan 26* (r) 



n G (r) + K e - 2fc 2 ' 



0; = refers to initial mixing angle at where neu- 
trinos are produced (i.e. at the neutrinosphere), and 
9f = 9{rf) is the mixing angle at rf where the neutrino 
helicity flip probabilities are evaluated. 

Finally, we note that the weak field approximation is 
fulfilled since 



AGM R 2 
5f 3 



>R< 1Q- 2 , 



as u> < 10 4 Hz, and that rotational effects are relevant 
during the time scale to < 10 sec (to is the time scale 
for the emission of the energy ~ 0.5 x 10 53 er g by each 
neutrinos degree of freedom with p ~ lOMeV) |33j . 



It has been suggested a mechanism for the generation 
of pulsar kicks which accounts for the magnetic and grav- 
itational fields of a rotating nascent star. Owing to the 
relative orientation of neutrino momenta with respect to 
the direction of the angular velocity, the energy flux turns 
out to be generated anisotropically. This effect is related 
to spin-gravity coupling, which differs for neutrinos with 
opposite helicity (left- and right-handed helicity). The 
estimations on Am 2 — sin 2 9 parameters agree with the 
best fit of the present data on neutrino oscillations. 

Some conclusive remarks are in order: 

• Results discussed in this contribution, as well as 
the papers \v\ fl8L IT^ | , in which pulsar kicks are 
discussed in relation to gravitational waves, have 
been obtained in scmiclassical approximation, i.e. 
the gravitational field is described by the classi- 
cal field equations of General Relativity. It will be 
of interest to investigate within the framework of 
quantum theories of gravity. 

• Results suggest a correlation between the motion 
of pulsars and their angular velocities. Such a con- 
nection seems to be corroborated by r ecent analysis 
and observations discussed in ji^, |5(| . 

• The mechanism here proposed is strictly related to 
the gravito-magnetic effect, an effect predicted by 
General Relativity [3(j, as well as by many metric 
theories [SJ. Its origin is due to the mass-energy 
currents (moving or rotating matter contribute to 
the gravitational fields, in analogy to the magnetic 
field of moving charges or magnetic dipole). Ex- 
periments involving the technology of laser ranged 
satellites [52| are at the moment the favorite can- 
didate to test gravitomagnetic effects. 

In connections with the mechanism proposed in this 
paper, a direct evidence of the gravitomagnetic ef- 
fect seems to be provided by pulsar kicks. Future 
investigations on the velocity distribution of pulsars 
will certainly allow to clarify this still open issue. 
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